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Proton acceleration by circularly polarized traveling electromagnetic wave 

Amol Holkundkar,* Gert Brodin, 1 " and Mattias Marklund* 

The acceleration of charged particles, producing collimated mono-energetic beams, over short distances holds 
the promise to offer new tools in medicine and diagnostics. Here, we consider a possible mechanism for accel- 
erating protons to high energies by using a phase-modulated circularly polarized electromagnetic wave propa- 
gating along a constant magnetic field. It is observed that a plane wave with dimensionless amplitude of 0.1 
is capable to accelerate a 1 KeV proton to 386 MeV under optimum conditions. Finally we discuss possible 
limitations of the acceleration scheme. 



I. INTRODUCTION 



II. MODEL DESCRIPTION 
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Laser induced particle acceleration has drawn considerable 
interest among researchers all over the world since the pio- 
neering work by Tajima and Dawson [1]. The acceleration 
gradient of conventional linear accelerators is of the order of 
10 5 V/cm, however today's state of the art lasers are capable to 
produce the acceleration gradient many orders of what can be 
achieved using conventional LINACs. In general laser based 
accelerators can be divided based on the medium in which 
the acceleration takes place, which can be either vacuum or 
a plasma. The vacuum as a medium for particle acceleration 
has some inherent advantages over plasma medium. The prob- 
lems like instabilities are absent in vacuum, it is easier to in- 
ject the pre-accelerated particle beam in vacuum as compared 
to the plasma, collisions of particles with media causing en- 
ergy loss and beam spreading is ruled out, etc. Thus we will 
focus on particle acceleration in vacuum in this article. 

The relativistic motion of the charge particle in large am- 
plitude electromagnetic (EM) fields are studied in detail by 
many authors. The motion of the charged particle in trans- 
verse EM wave and the constant magnetic field along wave 
propagation is studied by Roberts and Buchsbaum [2], which 
was further extended analytically and experimentally by Jorv 
and Trivelpiece [3]. More recently the in-depth Hamiltonian 
analysis of the dynamics of charge particle in a circularly po- 
larized traveling EM wave is been studied by Bourdier and 
Gond [4] . Various different schemes have been proposed for 
the acceleration of charge particle in traveling EM wave [5- 
8], some includes the homogeneous magnetic field however 
some includes the two counter propagating EM waves. 

In this paper we will consider an alternative method to ac- 
celerate protons in vacuum by circularly polarized electro- 
magnetic waves, where the main new ingredient is a phase- 
modulation of the EM wave Emphasis would be on under- 
standing the dynamics of proton motion under the proposed 
scheme. The next section will briefly describe the proposed 
scheme followed by the results and discussion. 
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A circularly polarized traveling wave propagating along the 
z direction is considered. The electric and magnetic fields of 
the wave are given by, 



E x = En sin [co (t -z/c) + (j>(t - z/c)} 



(1) 



E y = -Eq cos[co(t-z/c) + <p{t- z/c)] (2) 

and the magnetic fields are expressed as, 

B x = {E /c) cos[co(t-z/c) + <j>{t-z/c)] (3) 



By = (Eq/c) sin [co(t - z/c) + (f - z/c)} 



(4) 



with, <j>(t — z/c) is the phase modulation function which is 
given by 



<p(t-z/c) = % sin[aco(t—z/c)] 



(5) 



where, a is the so called phase modulation factor which con- 
trols the extent of the modulation. A constant magnetic field 
is also applied along the direction of wave propagation given 
by B z = b . 

The electric fields denoted by Eq. 1 and 2 can be generated 
by introducing an electro-optic phase modulator. This is an 
optical device in which a element displaying the electro-optic 
effect is used to modulate the beam of light. The modulation 
can be done in phase, frequency, polarization, and amplitude. 
The simplest kind of modulator consists of a crystal, such as 
Lithium niobate whose, refractive index is a function of the 
applied electric field [9]. An appropriate electric field along a 
crystal can be applied in such a way that its refractive index 
modulates, which eventually will introduce the phase lag in 
the laser beam. The phase modulation will also depend on 
the length of the crystal and other parameters. Designing an 
accurate phase modulator for a specific problem may be an 
engineering concern, but for the purpose of this article we will 
assume that such a problem can be solved satisfactorily. 

The schematic diagram for the proposed scheme is shown 
in Fig. 1. The laser pulse is initially passed through the phase 
modulator so that the spatial and temporal dependence of the 
electric and magnetic fields are modified according to Eqs. 1, 
2, 3 and 4. This modified pulse is then injected into the ac- 
celerating cavity, protons under the influence of this modified 
laser pulse undergoes the acceleration along the transverse di- 
rection. 



Phase modulated pulse 




Power Supply 
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FIG. 1. Schematic diagram of the proposed scheme. 



An exact analytical treatment of the problem seems to be 
too involved because of the nature of electric and magnetic 
field profiles. It would be a non-trivial task to solve the mo- 
mentum equations corresponding to the field equations (1)- 
(4). In view of this we have numerically analyzed the dy- 
namics of the particle under the influence of the given field 
profiles. 

As is wellknown, the motion of the relativistic particle is 
described by the following equations, 



p/mo 



v = 



y/l + |p|7('«oc) 2 



dv 
It 



(6) 

(7) 

(8) 



where, p, v, r and »?o are relativistic momentum, velocity, 
coordinate and mass of the particle. The above equations 
are solved numerically by a standard Boris leapfrog scheme 
where particle motion is decomposed into motion in the elec- 
tric field and the motion in the magnetic field [10]. The par- 
ticle orbits is calculated by substituting Eq. (1) -(4) into the 
equation pf motion, and specifying the initial condition for the 
injection energy, letting the initial velocity be directed along 
Z. In rest of the paper we have used the dimensionless units 
for all physical quantities. 

For all the results presented here, the amplitude of the cir- 
cularly polarized wave is considered to be an = 0.1 (unless 
otherwise stated), where «q = eE/m e coc [11]. Similarly mag- 
netic field is denoted by bo = eB/m e co. Here, e and m e being 
the charge and mass of electron, E and B is the amplitude 
of electric and magnetic field in SI units with ft) being the 
laser frequency. The dimensionless space and time are taken 
in units of the wave number k and the angular frequency ft) 
respectively. 



III. NUMERICAL ANALYSIS 

In this article our main focus is to understand the dynam- 
ics of a single proton with energy 1 KeV (unless otherwise 



stated), injected along the propagation direction (z) of a phase 
modulated circularly polarized wave with amplitude «o =0.1 
and a constant magnetic field with magnitude of bo — 1 .0 (Fig. 
1) along wave propagation. Although it will be clear later 
in the article that this scheme is equally valid for the proton 
beam having some energy spread. The central theme of the 
proposed scheme is the introduction of the so called phase 
modulation factor (a) which can be expressed as the ratio of 
the phase modulation frequency (©,„) to wave frequency ft) i.e. 
a = (Q m /(Q. Next we will see how the value of a affects the 
resulting dynamics of the particle. 

The time evolution of the transverse kinetic energy of the 
particle for the phase modulation factor a = l/n with n = 2, 3 
and 4 are presented in Fig. 2 (a). It is observed that the pro- 
posed scheme of acceleration works well only when the wave 
frequency is some harmonics of the modulation frequency. A 
large deviation from this condition destroy the acceleration 
mechanism completely, Fig. 2 (b). Furthermore, as can be 
seen from Fig. 2(a), the efficiency of the scheme deteriorate 
gradually with higher harmonics of the phase modulation fre- 
quency. 

It would be interesting to analyze the underlying principle 
for the success of this acceleration scheme. It turns out that 
the magnetic field as well as relativistic effects plays only a 
minor role in this regard. Thus, for the sake of simplicity let 
us consider the following equation of motion, 



dv x 

' dt 



qtE x 



(9) 



where m, and q is the mass and charge of the particle. Only the 
temporal dependence in Eq. 9 is considered for the purpose 
of analysis, which is valid assumption since the dynamics is 
mostly independent of the spatial coordinates. Now using the 
field profile given by Eq. 1 and phase modulation function 
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FIG. 2. Temporal evolution of transverse kinetic energy of proton for 
(a) three different phase modulation factor a = \jn with n = 2,3, and 
4 being an integer, similarly the case when n is not an integer is also 
presented (b). The amplitude of the EM wave an = 0. 1 and constant 
magnetic field along wave propagation bo = 1.0 is considered. 



(a) a = 0.25 



(b) a = 0.29 




FIG. 3. Numerical solution of Eq. 9 with different values of a (a) 
0.25, (b) 0.29, (c) 0.33 and (d) 0.42 is presented, considered. 



given by Eq. 5, this simple equation of motion (Eq. 9)can be 
integrated to give 



v ( t ) = ^El f sinW + nsin{ax')]dx' +C (10) 
'«; Jo 

where T is time in dimensionless units and C is a constant of 
integration. For simplification we have chosen qiE x /nii = 1 
and C = 0. The solution of the above integral (Eq. 10) for 
different values of a are shown in Fig. 3. It is quite clear that 
if a = l/«, with n being an integer the temporal evolution 
of velocity is very much regular and periodic (Fig. 3 a,c), 
however for the values when n is not an integer the solution of 
the integral is irregular, (Fig. 3 b,d). The observation is that 
the success and failure of the proposed acceleration scheme 
mostly depend on the property of the integral given by Eq. 10. 
It is observed that acceleration works only when the frequency 
of the electromagnetic wave is some harmonics of the phase 
modulation frequency, i.e. a = l/n with « being an integer. 

Now let us focus on the trajectories of the particle under the 
influence of the phase modulated EM wave (Eq. 1, 2, 3 and 4) 
and external magnetic field (B z = bo). The time evolution of 
the particle trajectory along 3 space dimensions is presented 
in Fig. 4 (a). Here we have considered the value of a = 0.5 
(the same results holds qualitatively for other valid values of 
a, i.e. a = 0.20, 0.25 and 0.33). 

The longitudinal displacement along z is mainly governed 
by the energy at which the particle is being injected in the 
cavity (magenta curve), however oscillatory motion along y 
direction (red curve) is because of the presence of the external 
magnetic field along z direction. 

As can be infered from the Fig. 4 (b) and (d), the biggest 
problem with the acceleration scheme is the excursion along 
the x direction which is much larger than the displacement 
along the z direction of propagation. In this scenario the 
proposed scheme would be impossible to implement in prac- 
tice because of the excursion along transverse direction is too 
large. In order to deal with this problem a small detuning pa- 



rameter 8 can be added to phase modulation factor a, which 
would be helpful to confine the particle orbits in the x-y plane. 

Let us examine the trajectory of the particle after the addi- 
tion of a small detuning parameter 8 = 10~ 4 which modifies 
the phase modulation factor a from 0.5000 to 0.5001. It can 
be observed from Fig. 5(a) that the kinetic energy of particle 
increased by a factor of about 2 as compared to case when no 
detuning parameter is present (Fig. 2). Furthermore the tra- 
jectory of the particle is modified significantly after the intro- 
duction of the detuning parameter. As can be seen, the particle 
trajectory in the x-y plane (Fig. 5(d)) is closed which make it 
possible for the acceleration scheme to work in practice. 

It should be noted that the selection of the detuning param- 
eter 8 is very crucial for the success or failure of the accel- 
eration mechanism. The value of 8 must be small enough to 
make the scheme work in favor of acceleration, i. e. we must 
still have the wave frequency to approximately be a harmonic 
of the phase modulation frequency, such that the foundation 
of the scheme is not destroyed, see Fig. 2 (b). Moreover, 8 
cannot be too small, as it must be large enough to prevent the 
large excursions (Fig. 4 b,d). 

The effect of the detuning parameter on the energetics of 
the particle is presented in Fig. 6. As can be seen from this 
figure, reducing the detuning parameter from 8 — 0.0008 to 
8 = 0.0001 increases the efficiency of the acceleration signif- 
icantly. On the other hand, the transverse kinetic energy of 
the particle is directly related to the transverse excursion of 
the particle. For lower values of the 8 the particle orbits are 
larger such that the electric field of the wave tend to do more 
work for each orbital motion of the particle. As the detuning 
increases the particle orbit becomes shorter and shorter, result- 
ing in lower gain in energy from the wave. Apparently there 
is an optimum detuning at which one is able to gain maximum 
energy per orbital cycle of the particle. 

So far we have presented all the results with constant mag- 
netic field of amplitude bo = 1 . It can be understood that the 






FIG. 4. Temporal evolution of the particle trajectory along three 
space dimensions (a) is plotted along with the trajectory in x-z plane 
(b), y-z plane (c) and x-y plane (d). Here, we have considered the 
phase modulation factor a = 0.5. 
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FIG. 5. Temporal evolution of the transverse kinetic energy (a) is 
plotted along with the trajectory in x-z plane (b), y-z plane (c) and 
x-y plane (d) with the detuned phase modulation factor a = 0.5001. 
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FIG. 7. Temporal evolution of the transverse kinetic energy for vari- 
ous values of applied magnetic fields with detuned phase modulation 
factor of 0.5001. 
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FIG. 6. Temporal evolution of the transverse kinetic energy for vari- 
ous values of detuning parameter 8. The resulting phase modulation 
factor would be a + 5 with a = 0.50. 



field bo = 0.3 is shown in Fig. 8. As can be seen from Fig. 8 
the excursion in the transverse direction is controlled. Thus in 
this case, particle acceleration using a phase modulated elec- 
tromagnetic wave propagating along a constant magnetic field 
seems to be possible in a real scenario. 



IV. DISCUSSION 

Here we have focussed on the understanding of the underly- 
ing dynamics of the acceleration. We have established that the 
phase modulated circularly polarized wave can accelerate par- 
ticles in the presence of a constant magnetic field. There are 
optimum conditions on the detuning parameter and the magni- 
tude of the magnetic field, in order to gain maximum energies. 
However, for the purpose of applications, the spectral proper- 
ties plays a crucial role, in addition to the maximum energy. 



magnetic field is also responsible along with the detuning pa- 
rameter to restrict the excursion of particle in transverse di- 
rection. In view of this it would be interesting to see how the 
acceleration efficiency varies with the applied constant mag- 
netic field. 

The temporal evolution of the energy varying the constant 
magnetic field is presented in Fig. 7. The detuned phase mod- 
ulation parameter in this case is chosen to be 0.5001. It is 
observed that the proton can now be accelerated to energies 
of about 386 MeV when the applied magnetic field strength 
is 0.30. Furthermore, there seems to be an optimum magnetic 
field for the maximum acceleration of the particle. This be- 
havior can be explained on the basis of the Larmor radius of 
the particle which varies with the applied magnetic field. The 
magnetic field should be strong enough to bend the particle 
to avoid large excursion and should not be so large such that 
the particle orbit is very small and the resulting energy gain in 
one orbital motion is small. The particle orbit for the magnetic 




FIG. 8. Particle trajectory under the influence of applied magnetic 
field of magnitude 0.3. The particle is propagating along the z direc- 
tion. 
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FIG. 9. Temporal evolution of the kinetic energy for different values 
of injection energies (a). A energy spectrum of the output beam is 
calculated by injecting a beam with Gaussian energy spectrum (b). 
The inset in (b) represents the output beam energy as a function 
of injection energy, circles denotes the numerically measured val- 
ues however solid line is 10 order polynomial fit which is used in 
calculating the energy distribution function of output beam. 



In view of this it would be important to understand how this 
scheme works for a particle beam with an energy spread, in- 
stead of just a single particle. In order to shed some light on 
the acceleration dynamics of the particle beam we have sepa- 
rately simulated the single particle motion when injected with 
different energies, which can be perceived as the energy spec- 
tral spread of the beam. 

The time evolution of the transverse kinetic energy of sin- 
gle particle is presented in Fig. 9(a) for different injection 
energies. The phase modulation factor is considered to be 
a = 0.5001 and the magnitude 0.5 is considered for the con- 
stant magnetic field. It can be observed that even if the in- 
jection energies are varied an order of magnitude, the output 
energy of the beam is not very much affected, keeping the 
scheme functional. The most important property of the parti- 
cle beam is its energy distribution function. The energy distri- 
bution function of the output beam is calculated in Fig. 9(b) 
by considering the input beam having Gaussian energy spec- 
trum with peak energy of 600 KeV and FWHM of 166 KeV, 



the energy spread (AE/Ep ea k) of the input beam is about 28%. 
The peak energy of the output beam is observed to be 167 
MeV with FWHM of 5 MeV, with the energy spread about 
3%. 

The inset in Fig. 9(b) shows the variation of output energy 
with the injected particle energy. Circles denotes the actual 
numerical value which is fitted with a 10' /? order polynomial 
(solid line) in order to find the energy spectrum. The decline in 
the output particle energy with input energy can be explained 
on the basis of the interaction time of the particle with the 
EM wave. The faster the particle, the lesser is the interaction 
time with the EM wave, and hence the energy transfer to the 
particle. 

As we have observed in Fig. 9(b), the output particle en- 
ergy is more or less independent of the injected particle en- 
ergy, which is apparently visible in the almost monoenergetic 
energy distribution of the output beam. 



V. FINAL REMARKS 

A further evaluation of the feasibility of the proposed ac- 
celeration scheme need to go beyond the ID- variations of the 
electromagnetic fields. In particular, it is clear that picking pa- 
rameter values corresponding to intense lasers give succesful 
results for the maximum particle energies. However, a poten- 
tial limitation is that large transverse particle excursions ex- 
clude the use of very focused pulses, in which case the avail- 
able laser intensity drops accordingly. To some extent this 
may be remedied with a strong value of the static magnetic 
field, which limits the transverse excursion. In practice, mag- 
netic field stengths well beyond 100 T is needed if the system 
works in the optical laser regime, which makes this regime 
less attractive. Decreasing the wave frequency in the scheme 
reduces the need for extreme magnetic field strengths, since 
we may allow for somewhat larger particle excursions. The 
optimal frequency regime may lie in the infra-red regime or 
lower, but a full 3D analysis is needed to optimize the param- 
eters in a realistic scenario. This remains a project for further 
research. 
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